www.daoe-jotirnal.org 


Development and Applications of Oceanic Engineering (DAOE) Volume 3, 2014 


Thermohaline Convection in a Porous 
Medium in the Presence of Magnetic Field 
and Rotation 

Abdullah A. Abdullah, Sahar Z. Alkazmi 

Department of Mathematical Sciences, Umm Al-Qura University, P. O. Box 6337, Makkah, Saudi Arabia 
aamohammad@uqu.edu.sa 

Received 26 January, 2014; Revised 10 March, 2014; Accepted 20 March, 2014; Published 18 May, 2014 
© 2014 Science and Engineering Publishing Company 


Abstract 

The onset of Benard convection in a horizontal porous layer 
permeated by an incompressible, thermally and electrically 
conducting solute fluid under the effect of a uniform vertical 
magnetic field and a uniform vertical rotation is investigated. 
The porous layer is assumed to be governed by the Brinkman 
model. Analytical and numerical solutions are obtained for 
the cases of stationary convection and overstability. The 
critical thermal Rayleigh numbers are obtained for different 
values of the permeability of porous medium, 
Chandrasekhar number and Taylor number for different 
cases of solute boundary conditions. The related eigenvalue 
problem is solved using the Chebyshev polynomial Tau 
method. 
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Introduction 

Thermal instability theory has attracted considerable 
interest and has been recognized as a problem of 
fundamental importance in many fields of fluid 
dynamics. This theory has been enlarged by the 
interest in hydrodynamic flows of electrically 
conducting fluids in the presence of magnetic field and 
rotation (Chandrasekhar 1960). 

The onset of convection induced by thermal and solute 
concentration gradients in a horizontal layer of a 
viscous fluid is studied by Nield (1967) when the layer 
is heated from below and the solute concentration 
increases upwards for general set of boundary 
conditions. Nield (1968) studied the thermohaline 
convection in a porous medium for a general set of 
boundary conditions when the fluid is heated and 


soluted from below. The thermosolutal convection in a 
fluid layers heated and soluted from below in a porous 
medium has been studied by Sharma and Kumari 
(1992) to include, separately, the effects of a uniform 
horizontal magnetic field and a uniform rotation when 
the boundaries are free. Khare and Sahai (1992) 
studied the thermosolutal convection in a 
heterogeneous fluid layer heated and soluted from 
below in porous medium when the boundaries are free 
using Brinkman - Oberbeck - Boussinesq equations for 
the cases of stationary and overstability convection. 
This problem has been extended by Khare and Sahai 
(1993) to include the effect of magnetic field and by 
Khare and Sahai (1994) to include the combined effects 
of rotation and magnetic field. The thermosolutal 
convection of a horizontal layer of an incompressible 
viscous fluid in the presence of vertical magnetic field 
is discussed by Al-Aidrous and Abdullah (2001). 
Related problems on the thermosolutal convection in a 
porous medium have been discussed by Sunil et al. 
(2005), Singh and Kumar (2011), Nield and Bejan (2013) 
and others. 

This work studies convective instability in a horizontal 
porous layer permeated by an incompressible, 
thermally and electrically conducting solute fluid 
using Brinkman model in the presence of a uniform 
vertical magnetic field and a uniform vertical rotation. 
Analytical solutions were obtained when both 
boundaries are free and numerical results were 
presented for the cases of free and rigid boundaries. 
The numerical computations were performed using 
the Chebyshev Tau method. 

The implementation of Chebyshev Tau method is 
highly useful in obtaining accurate eigenvalues for one 
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layer and two layers problems. The application of a 
D~ Chebyshev tau methods to a variety of subject 
areas is discussed by Abdullah and Lindsay (1990, 
1991a, 1991b), Lindsay and Ogden (1992), Abdullah 
(2000), Straughan (2001, 2002), Al-Aidrous and 
Abdullah (2005), Gheorghiu (2007), Gheorghiu and 
Dragomirescu (2009), Banjar and Abdullah (2010) and 
others. 


Mathematical Formulation 


Consider an infinite horizontal layer of porous 
medium confined between two parallel horizontal 
boundaries and permeated by an incompressible 
thermally and electrically conducting viscous solute 
fluid. The porous layer is subjected to a constant 
vertical magnetic field with a magnetic field intensity 
B (=(0,0,B)) and to a rotation about the x 3 axis with 
angular velocity Q . Gravity g acts in the negative 
direction of x 3 . 

Let 6 be the absolute temperature of the fluid C the 
solute mass concentration, T 0 , T 0 -T , T > 0 , the 
temperatures at x 3 =0 and x 3 =d respectively, 
S 0 , S 0 -S the concentration at x 3 = 0 and x 3 =d 
respectively, p t , the density of the fluid at T 0 and a , 
a' the coefficients of volume and solute expansion 
respectively, then the fluid density is proportionate to 
0 and C such that 

p = p 0 (\-a9 + a'C). 


Thus the governing equations are 

v i = -( P / p 0 ) i + vV 2 p -g(l-a9 + cc'C)S i3 + —B ik B k 

PoB 

~^~ u i +2 e ijk v P k > 

°i,i = 0 , 

9 = K V 2 # , 

C = k'S/ 2 C , 


Bu= 0 , 


'f ! ^ij k H k ,j 


dB, 

dt 


L = o iJ -VjB iJ -j 1 e i j k J kj 


(1) 


where p, H i and P are velocity, magnetic field and 
pressure, v, p, k l , k, k and rj are kinematic viscosity, 
magnetic permeability, permeability of porous 


medium, coefficient of thermal diffusivity, coefficient 
of solute diffusivity and electrical resistivity. Equations 
(1) have a steady state solution of the form 

p = 0 , 7,. =0,0 = 9{x 3 ) = T a -p x 3 , P = P(x 3 ) 

C = C (x 3 ) = S a - P'x 3 , 5, = (0, 0, B) B = constant (2) 


where y 6=T/d is the uniform adverse temperature 
gradient and / 8'=S/d is the uniform adverse 
concentration gradient. Suppose now that the variables 
in (2) have a small perturbation of the form 

p = 0 + s p, 7, =0 + £ 7 j,0 = 0 + s 9 ,C = C + s C , 

P = P + £ P, B; =(0,0 ,B) + eb, 

where p , 7, , 0, C, P, b j are the linear perturbations of 
velocity, current density, temperature, solute 
concentration, pressure and magnetic induction 
respectively. Following standard procedures the linear 
equations of (1) become 

~rr = ~Pj + V 2 p + JR, 98, 3 - JJ CS i3 ~p + b u3 
dt N 

+ Vr e ijk Uj S k3 
P,=0 , 

p — + h Sr u 3 =v 2 9 , 
r dt V ' 3 

p' r ^L + H'jR s p=V 2 C , (3) 

dt 

hi - o , 
h ~ e ijk b k j 

db , 

h - Q°i,3 - e ijk J k,j • 

Note that 


R, = 


d \ R ,L*1 


K V 


K V 


d\ 



P V 

n 


<2 = 


B 2 d 2 
Po V ’ 




T = 


4 n 2 d 2 


v 

K 


where R , , R, are the thermal and solute Rayleigh 
numbers, P r , P m are the viscous and magnetic 
Prandtle numbers, P m is Schmidt number, Q is the 
Chandrasekhar number, T is the Taylor number and 

| 1 when heating from above 
[-1 when heating from below 
1 when solute increases above 
-1 when solute increases below 


H = 


H' = r/r 


\P'\ 

1 1 
F 
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Boundary Conditions 

The fluid layer is confined between the planes x 3 = 0 
and x 3 = d and on these planes we need to specify 
mechanical, thermal and magnetic conditions. 
Mechanical conditions mean rigid or free boundaries. 
For rigid boundaries 

u 33 =0, = 0 on x 3 = 0, d , 

and for free boundaries 

u 3 , 33 = 0 ’ £ 3,3 0n x 3 =0,d. 

where c is third component of the vorticity. Thermal 
conditions mean an insulating or a conducting 
boundary. For a conducting boundary 6=0 ext and for 

an insulating boundary x =0, where 0 ext is the 

temperature outside the conducting boundary. 
Electromagnetic conditions assume either a stationary 
perfectly conducting electromagnetic boundary in 
which b 3 = 0 or a perfectly insulating electromagnetic 
boundary in which (curl B) 3 = 0 . Solute conditions 
assume either an impermeable boundary in which 
C 3 =0 or a permeable boundary in which C=C ext 

where C ext is the solute concentration in the exterior 


L=D 2 -a 2 . We may eliminate J,b,6,%,C from these 
equations to obtain a 14* order ordinary differential 
equation in iv of the form 

L(L- a)(L - aP )(L - aP r )(L - aP' r )[(L - aP m )(L -a~) 

N 

- QD 2 ]w-a 2 HR,(L- aP m )(L - aP,!)[(L - aP m )(L — u — 

N 

-QD 2 ]w + a 2 H' R (L- a P m )(L - aP r )[(L - aP m )(L - a - -J-) 

N 

-QD 2 ]w- QL(L -aP)(L- aP' r ) [(L - aP m )(L - cj - -J-) 

N 

-QD 2 ]D 2 w - -J- L(L - aP m )(L - aP r )(L - aP;.)[(L 
N 

-aP m )(L-cj-±-)-QD 2 ]w+ T(L - aP ) 2 (L- crP )(L 
N 

-aP;.)D 2 w = 0. 

( 5 ) 


The Free Boundary Problem 

Suppose that both boundaries are free and 
w = A sin( 1 7i x 3 ), where A is a constant and 1 is an integer. 

Thus Lw = -Aw ; A = l 2 n 2 + a 2 and equation (5) 
becomes a sixth order polynomial in a with real 
coefficients. Now we have four cases to consider. 


region. 


The Eigenvalue Problem 

Now we shall construct the related eigenvalue 
problem from equations (3) and the boundary 
conditions. We now apply a normal mode expansion of 
the form 


(p = (pixQex-pidnx^ +mx 2 ) + at \ 


where n,m are the wave numbers of the harmonic 
disturbance and a is the growth rate. Thus equations 
(4) become 


<r E = Ld + DJ - — <f + Vr Dw, 

N 

a Lw = l}w-a 2 jR, d + a 2 x Jr^ C + UDb) 

- — Lw - ~Jt Dd, 

N 

aP m J = QD% + LJ, (4) 

a P m b = QDw + Lb, 
a P r e = Le-H^R, w, 
a P X C = LC-H' JR S w. 

where a = yj n 2 +m 2 is the wave number, D = d/C and 


Case 1. The porous layer is heated from above and the 
solute concentration decreases upwards 

Here H = 1, //' = -! . For stationary convection case we 
set er = 0 in equation (5) to obtain the thermal 
Rayleigh number 

R,=-[-4- (U 2 +Ttt 2 A) + R s ] 
a~C 

A 

where U = (A 2 + — + Q n 1 ) from which we deduce 
N 

that no stationary convection stability. To study the 
possibility of overstability suppose that ^ 0 and we 
obtain from equation (5) 

D y (y + iaQiy + iaf,) (r + i^P,.) 

K\ = — A, 

O-l) (r + P r ) 

r (y + i<j\P r ) y (r+icTjP,.) j_ 

(r-m+i^pj 1 0-1) N, 

, P+'<h^„)P+'<hQ) T 

ir - 0[(/ + icr ] P m )(y + fcj| + ' ) - 2, ] 

h'l 

where y = 1 + a, . To study the effect of solute Rayleigh 
number, magnetic field, permeability of porous 
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medium and rotation on the thermal Rayleigh number 

, , . , dR n dR n dR„ dR n 

we have to discuss the nature of — — , — — , — — , — — 

dR sl dQ, dN, dT, 

In each of these derivatives we separate the real and 
imaginary parts, then solve the resulting equations to 
obtain : 

dR t dR t y dR t y 2 

dR s ' dQ, (J-I)' dN, N 2 (y—l ) ' 

dR, _ y 2 +a , 2 P m P r 

dT i {y-l)[y{y+HN,)+a, 2 P m -Q,] 

Thus all the parameters have a destabilizing effect. In 
fact rotation could have a stabilizing effect provided 

Qi>r(r+^r)+^ p ,n 

Case 2. The porous layer is heated from above and 
the solute concentration increases upivards 

Here H = 1, H' = 1. For stationary convection case we 
set cr = 0 in equation (5) to obtain 


From which we obtain 

dR n dR n y dR,, y 2 

dR sl dQ, a, ’ dN, a, N 2 ’ 

dR, i y 2 +a?P m P r 

dT, 1, 

a \ r(r+Y ~Qi 

Thus solute Rayleigh number and permeability of 
porous medium have a stabilizing effect whereas 
magnetic field has a destabilizing effect. Moreover 
rotation has a stabilizing effect provided 

Case 3. The porous layer is heated from below and the 
solute concentration increases upivards 

Here H =— 1, H' = \ . For stationary convection case we 
set a = 0 in equation (5) to obtain 

R, = —j— (U z +Ttt 2 A)-R s . 

a 2 U 


R, (U 2 + Tx 2 A)] + R s . (6) 

a~U 

Note that stationary convection stability is possible 

provided R s >-^ — {U 2 + Tk 2 A). 
a~U 

To study the effect of solute Rayleigh number, rotation, 
magnetic field and permeability of porous medium on 
R t we deduce from (6) that 


dR, _ ! dR, _ _ TV dR, _ _ n 2 A U 2 -tt 4 A 2 T 
dR s ’ dT a 2 U ’ dQ a 2 U 2 

\\{-U 2 +Tk 2 A) 
dR, _ jy 2 

dN ~ a 2 U 2 

Thus the solute Rayleigh number has a stabilizing 
effect, whereas rotation has a destabilizing effect. 
Magnetic field and permeability of porous medium 
have a destabilizing effect provided An 2 T>U 2 . For 
overstability case we have 


y (y + icr, )(y + ia,P r ) | (y + icr,P r ) 

(y-l) (y + icr,P ' r ) !l 

y ( y + icr x P r ) Q r (y + ia,P r ) J_ 

(y-Y)(y + icr l P m ) 1 (y-l) N, 

(r + iviPJir + i°\P r ) 

l 1 

{y-\)[{y+KT l P m ){y+i(J l +— )-£>,] 


Clearly stationary convection stability is possible 
2 

whenever R <— — (U 2 +Tx 2 A) and 
a 2 U 

dR, _ dR, _ A 2 k 2 dR, _ ji 2 AU 2 - n A A 2 T 
dR s ’ dT a 2 U ’ dQ a 2 U 2 

/p ^(-U 2 +tt 2 AT) 
dR, N 2 

dN crU 2 

Thus solute Rayleigh number has a destabilizing effect, 
whereas rotation has a stabilizing effect. Magnetic field 
and permeability of porous medium have a 
destabilizing effect provided An 2 T>U 2 and 

permeability of porous medium has a stabilizing effect 
provided An 2 T >U 2 . For overstability case 

R _ 7 (r + io-iXr + ioV) (y + ia,P r ) R 
n (y-l) (y + i<j,P,') sl 

| y (y + ia,P r ) | y (y + ia,P r ) J_ 

(y-\)(y + i<j,P m ) ^ (y-l) N, 

(r + 'oV,)(r + 'oV) T 

(r - i)l(y + )(y + *./)' - Q, 1 

From which we obtain 

dR n = _j dffi = 7_ fffn = 7' 

dR sl ’ dQ, a , dN , a, N 2 ’ 

dRn 7 2 +cT{P m P r 

dTl r(r+^)+<rfP m - Qi 
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i.e. solute Rayleigh number and permeability of porous 
medium have a destabilizing effect, whereas magnetic 
field has a stabilizing effect. Rotation has a stabilizing 
effect provided 

Qi>rtr + -^-) + °i P^- 

Case 4. The porous layer is heated from below and the 
solute concentration decreases upivards 

Here H = — 1, H' = — 1. For stationary convection case 
we set cr = 0 in equation (5) to obtain 

R, = (U 2 +T Ti 2 X) + R 

a 2 U 

and hence 

dR, _ l dR, _ X 2 n 2 dR, _ tcX U 2 -n 4 X 2 T 
dR s ’ dT a 2 U ’ dQ a 2 U 2 

, R ^(-U 2 +x 2 AT) 
dR t _ N 2 

dN a 2 U 2 

Thus the solute Rayleigh number and rotation have a 
stabilizing effect on the system. However magnetic 
field has a stabilizing effect when U 2 >A7CT and the 
permeability of porous medium has a stabilizing effect 
when A n 2 T > U 2 . For the overstability case 

D r (y + icrfiy + ia^f) ( y + ia \P r ) 

K,. = H n rl 

(r- 1) (r + i&t K ) 

| y(y + ia { P r ) | y (y + icr { P r ) J_ 

(r-m+i^PJ 1 (r-1) N t 

(r + iviP„,)(r+iviP r ) T 

(r~])[(r + ia i p m )(r+ia l + ' )-g,] 

from which we obtain 

= 1 = dR n _ y 2 

dR sl dQ l a, ' dN t aj N 2 

dR a r +<J{P,„P r 

dT, r 1 , 

a i r(r+-)+<j{P m -Q 1 

”1 

Thus the solute Rayleigh number and the magnetic 
field have a stabilizing effect whereas the permeability 
of porous medium has a destabilizing effect on the 
system. Rotation has a stabilizing effect on the system 
provided 

Q > r(.r+-^-) + v;p m . 


Numerical Results and Discussion 

The system of equations (5) together with the 
boundary conditions constitute an eigenvalue problem 
of order 14. This eigenvalue problem is solved using 
the Chebyshev Tau method when the fluid layer is 
heated from below and the solute concentration 
increases upwards which corresponds to case (3) in the 
previous section. 

Numerical solutions are obtained for free and rigid 
boundaries for the following cases 

(i) DC = 0 at X 3 = 0 , 1 l 

(ii) C= 0 at xt ) =Q and DC = 0 at X 3 =D 

Stationary Convection 

The relation between the magnetic parameter 
(Chandrasekhar number), Q, and the critical thermal 
Rayleigh number, Rt, is displayed in Fig. (1) when both 
boundaries are free and 7’ = 10 5 . As Q increases Rt 
increases for all values of N which indicates that 
magnetic field has a stabilizing effect. Moreover the 
values of Rt for case (i) of solute boundary conditions 
are less than the corresponding values of case (ii) of the 
conditions. In fact the difference between these values 
is very small when N=0.1 but when N=0.001 this 
difference is significant. 

Figs. (1) shows also that as the permeability of the 
porous medium, N, decreases Rt increases which 
indicates that N has a destabilizing effect on the system. 
Moreover numerical results show that as T increases Rt 
increases for all values of Q and N which indicates that 
rotation has a stabilizing effect on the system. Similar 
effects are obtained when both boundaries are rigid. 
These effects are displayed in Fig. (2). 

Overstability Case 

The relation between the magnetic parameter, Q, and 
the critical thermal Rayleigh number, Rt , is displayed 
in Figs. (3) when both boundaries are free and T = 10 5 
respectively. As Q increases Rt increases for all values 
of T when N=0.001 which indicates that magnetic field 
has a stabilizing effect. However as the permeability of 
porous medium increases unexpected results are 
obtained. In fact for small values of Q, as Q increases Rt 
decreases but for large values of Q as Q increases Rt 
increases as Q increases. These results coincide with 
those of Chandrasekhar (1981) where he considered 
the same problem in absence of porous medium. 

Moreover for small values of T the values of Rt for case 
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(i) of solute boundary conditions are greater than the 
corresponding values of case (ii) of the conditions 
when N=0.1, however when N=0.001 the values of Rt 
for case (i) of solute boundary conditions are less. But 
for large values of T this conclusion is reversed. 

Figs (3) shows also that as the permeability of the 
porous medium, N, decreases Rt increases which 
indicates that N has a destabilizing effect on the system. 
Moreover as T increases Rt increases for all values of Q 
and N which indicates that T has a stabilizing effect on 
the system. 

We also note that the case of overstability do not 
appear unless Pm> Pr and Q exceeds a certain critical 
value. This critical value is calculated when N = 0.1 and 
it is approximately equal to 1606 when 2=0, 1568 when 
T= 10 3 , 1366 when T = 10 4 and 380 when T= to 6 . It 
appears that this critical value decreases as T increases. 
Similar effects are obtained when both boundaries are 
rigid. These effects are displayed in Fig. (4). 

Fig. (5) shows the critical values of Q for T=50000 when 
N = 0.1 and N=0.001 for case (1) of solute boundary 
conditions when both boundaries are free.. The critical 
value of Q at which overstability is attained is 
calculated when N=0.001. This value is approximately 
equal to 2305 when 7=0, 2350 when T=50000 and 2400 
when T = 10 5 . 

Conclusion 

The thermohaline convection in a horizontal porous 
layer permeated by an incompressible solute fluid in 
the presence of magnetic field and rotation has been 
examined. Analytical solutions are obtained for 
different cases when both boundaries are free. 
Numerical solutions are obtained when both 
boundaries are free and rigid for the cases of stationary 
convection and overstability when the porous layer is 
heated from below and solute concentration increases 
upwards. 

For the stationary convection and overstability cases, 
the presence of both magnetic field and rotation 
stabilize the system for different values of the 
permeability of porous medium. It appears also that as 
the fluid becomes less porous the fluid becomes more 
stable. 
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